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REMARKS 

Claims 1-16 are now present in this application, with new 
claims 12-16 being added by the present Preliminary Amendment. 
It should be noted that the amendments to original claims 1-11 
of the present application are non-narrowing amendments, made 
solely to place the claims in proper form for U.S. practice 
and not to overcome any prior art or for any other statutory 
considerations. For example, amendments have been made to 
broaden the claims; remove reference numerals in the claims; 
remove/change any phrases unique to European practice; remove 
multiple dependencies in the claims; and to place claims in a 
more recognizable U.S. form, including the use of the 
transitional phrase "comprising" as well as the phrase 
"wherein" . Other such non-narrowing amendments include 

placing apparatus- type claims (setting forth elements in 
separate paragraphs) and method- type claims (setting forth 
elements in separate paragraphs) in a more recognizable U.S. 
form. Again, all amendments are non-narrowing and have been 
made solely to place the claims in proper form for U.S. 
practice and not to overcome any prior art or for any other 
statutory considerations . 
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CONCLUSION 

Accordingly, in view of the above amendments and remarks, 
an early indication of the allowability of each of claims 1-16 
in connection with the present application is earnestly 
solicited. 

Should there be any outstanding matters that need to be 
resolved in the present application, the Examiner is 
respectfully requested to contact Donald J. Daley at the 
telephone number of the undersigned below. 

If necessary, the Commissioner is hereby authorized in 
this, concurrent, and future replies, to charge payment or 
credit any overpayment to Deposit Account No. 08-0750 for any 
additional fees required under 37 C.F.R. § 1.16 or under 37 
C.F.R. § 1.17; particularly, extension of time fees. 



Respectfully submitted, 

HARNESS, DICKEY & PIERCE, P.L.C. 




P.O. Box 8910 
Reston, Virginia 20195 
(703) 668-8000 

DJD/bof 
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SUBSTITUTE SPECIFICATION 

Description 

VOXEL-DRIVEN SPIRAL RECONSTRUCTION FOR CONE BEAM COMPUTED 

TOMOGRAPHY 

[0001] This application is the national phase under 3 5 

U.S.C. § 371 of PCT International Application No. 
PCT/EP2004/000754 which has an International filing date of 
January 28, 2 0 04, which designated the United States of America 
and which claims priority on German Patent Application number 
DE 103 04 662.3 filed February 5, 2003, the entire contents of 
which are hereby incorporated herein by reference. 



Field 



[0002] The invention generally relates to a method for 
generating images in computed tomography with the aid of 3D 
image reconstruction metho d. For example, it relates to one— 
wherein^ in order to scan an examination object with the aid of 
a conical beam emanating from a focus and the aid of a planar, 
preferably multirow, detector for detecting the beam, the focus 
is moved on a spiral focal track about the examination object— 
. T- ferhe detector aupplying supplies output data, that correspond 
to the detected radiation, and image voxels from the scanned 
examination object are reconstructed from the possibly 
preprocessed output data and reproduce the attenuation 
coefficients of the respective voxel. 

Background 

[0003] Beams having cone beam geometry are uocd — ±& — fefee — prior 
art — a #known for computed tomography with multirow detectors. 
Reconstructing an image volume requires consideration of the 
cone beam geometry in the 3D image reconstruction, the cone 
beam geometry substantially complicating the reconstruction 
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problem. Two different groups of image generation methods are 
known in principle, specifically approximate methods and exact 
methods : 

Approximate methods : 

[0004] Approximate methods are distinguished by a high measure 
of practicability and flexibility. The angle of inclination of 
the measuring beams to the axis of rotation (cone angle) is 
considered in an approximate fashion, for which reason the 
error owing to the approximation grows with the cone angle. 
Starting from a certain number of detector rows, each 
approximate method will cause — image artof acto artif acts . A 
distinction is made between 2D and 3D methods in the case of 
the approximate methods . 

[0005] ASSR [1], AMPR [2] and SMPR [3] may be named as examples 
of the 2D method or 2D rebinning. In this case, synthetic 
projection data are approximated from the cone beam data, the 
geometry of the synthetic projection data being selected such 
that all the synthetic measuring beams lie in a flat plane. The 
object distribution in the flat plane can be reconstructed from 
the synthetic pro j ection data with the aid of the conventional 
algorithms of 2D CT . However, with growing cone angles the 
approximation in the rebinning step leads relatively quickly to 
image artof acto artif acts . 

[0006] In the case of the 3D method, filtering of the 
projection data and a subsequent 3D back projection are 
undertaken. The recorded geometry of cone beam type is taken 
into consideration here exactly in the back projection. Various 
feasible approximations are applied in the filter step. These 
approximation methods set forth below have in common the fact 
that the filtering consists of a ID ramp filter that is aligned 
in the direction of the projection of the spiral tangent. The 
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individual approximation methods differ in the processing of 
redundant data. 

[0007] On the one hand, it is possible to average in the axial 
direction, as is known from [4] . Here, after the filter step 
the filtered projection data are back projected into the 
reconstruction volume such that the measured data are 
accumulated in weighted fashion in the process. The weights 
result from the axial distance of the measuring beams from the 
voxel to be reconstructed. 

[0008] On the other hand, a reconstruction of axial slices is 
known from [5] . For this purpose, all the measured data that 

cut the image slice are used for the reconstruction. Data 

redundancies are considered approximately, neglecting the_— cone 
angle . 

[0009] Use may be made for this purpose of conventional methods 
such as, for example, the Parker weighting known from [6] . The 
filter step conoioto in this case of includes a weighting of the 
projection data in accordance with the present data redundancy, 
followed by convolution using a ID ramp filter. This method 
requires a relatively large detector surface, since the 
projection of the entire reconstruction slice must be included 
on the detector for each focal position. 

Exact methods: 

[0010] A comparative study of the most important exact 
algorithms is to be found in [7] . Exact methods consider the 
recorded geometry of cone beam type in a fashion free from 
error both in the filter step and in the 3D back projection. 
These methods achieve good image results independently of the 
cone angle that occurs. However, they are extremely complicated 
and very inflexible in application. For example, a reduced 
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pitch accompanied by the employment of data redundancies can be 
realized only to a limited extent, and there is no possibility 
at all of selecting measured data for cardiac imaging with a 
high time resolution. 

[0011] The problem now conaiato — resides in that the image 
quality achieved and the required computational outlay are 
oppositely oriented as a function of the method used. 
Algorithms having excellent image quality in conjunction with 
large cone angles, that is to say for detectors with many 
detector rows, are essentially intended for future scanner 
generations. It would be desirable to avoid exact methods, 
since these are complicated and inflexible. 

SUMMARY 

[0012] It is therefore — an object of an embodiment of the 
invention to find in computed tomography^ having multirow 
detectors and beams with cone beam geometry^ a proposal to 
improve or even solve the generation of images that, on the one 
hand constitutes an improved or even optimum compromise between 
image quality and computational outlay, and on the other hand 
can also be applied sufficiently flexibly. 

[0013] — Thia — obj cot — is — achieved by mcana — e€ — fehe — f caturco — e£ — fehe 
independent — patent — claimo . — Advantageous — dove 1 opmcn t o — eS — fefee 
invention arc the aubject matter of dependent claimo . 

[0013] The inventors have recognized the following: 

[0014] Starting from a reconstruction method such as is 
described in [5] , instead of reconstructing entire image slices 
it is also possible to reconstruct the image voxels separately. 

-tfe a-s feThen^ neceooary i** thio caoc fee proccoa data 

redundancies are processed by using methods that permit every 
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arbitrary projection angular range that is greater than the 
required half revolution n. 

[0015] What aspect which is novo 1 — unique in at least one 
embodiment of the invention is the separate reconstruction of 
individual voxels. The projection data are weighted separately 
for approximate consideration of data redundancies for each 
individual voxel. The result of this is that, by contrast with 
the conventional back projection, instead of a filtered 
measured value being smeared along a line over all the voxels 
cut by the line, the measured value filtered for the voxel V is 
accumulated only on the voxel V. 

[0016] i Fhe — At least one embodiment of the p resent invention 
offers at least one of the following advantages-? — : good image 
quality is expected even for large cone angles, flexibility in 
pitch up to the theoretical maximum value, use of all the data 
redundancies and the possibility of data selection, for 
example, for cardiac imaging with high time resolution. 

[0017] In accordance with this basic idea, in an embodiment of 
the invention the inventors propose the following method for 
generating images in computed tomography with the aid of a 3D 
image reconstruction method that has at least the following 
method steps: 

in order to scan an examination object with the aid of a 
conical beam emanating from a focus and the aid of a 
planar, preferably multirow, detector for detecting the 
beam, the focus is moved on a spiral focal track about the 
examination object, the detector supplying output data, 
that correspond to the detected radiation, and 
image voxels from the scanned examination object are 
reconstructed from the possibly preprocessed output data 
and reproduce the attenuation coefficients of the 
respective voxel, 
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each image voxel being reconstructed separately from 
projection data that comprise a projection angular range 
of at least 180°, and 

an approximate weighting taking place for each voxel 
considered in order to normalize the projection data used 
relating to the voxel. 

[0018] It is advantageous in at least one embodiment, 4r& — thio 
CQ3C — when - as illustrated in Figure 4 - in order to 
reconstruct an image voxel V, use is made of all the detector 
data along a straight line that runs through the cone beam 
projection of the voxel V and is aligned in the direction of 
the projection t of the spiral tangent S t . An improved or even 
optimum detector employment is achieved thereby. 

[ 0019 ] Furthermore , the image data of the detector image can be 
subjected to cosine weighting in order thereby to achieve a 
compensation of the oblicjue irradiation of the focus onto the 
detector . 

[0020] If the scanning operation does not have a sufficiently 
high resolution in order for all the required detector data to 
be obtained directly for a voxel being considered, data not 
directly available can be obtained from the available data by 
interpolation from neighboring detector data (detector pixel 
values) . 

[ 0021 ] Furthermore , during the weighting for compensating the 
data redundancy (normalization) two measuring beams (Sa, Sb) 
can be regarded as redundant precisely when it holds for the 
parallel coordinates 9, p that: (0 a = 2k-7i + 0 b and p a = Pt>) or 

(9 a = (2k + l)-7t + 9 b and p a = -p b ) , where k represents an 
arbitrary natural number, and 9 represents the projection 
angle, and the projection angle 9 represents the sum of focal 
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angle a and fan angle |3, and p represents the distance from the 
z axis. 

[0022] The redundantly available data can be multiplied by 
generalized Parker weights for the purpose of normalization. 
The normalized data are processed with the aid of a ramp filter 
modified by an advantageous smoothing window. It is 
advantageous, furthermore, when a distance weighting is 
performed for the 3D back projection into the voxel considered. 

[0023] Finally, it is to be pointed out that the method 
according to at least one embodiment of the invention can also 
be used for cardiac computed tomography by selecting, weighting 
or sorting in a manner known per se the measured data in 
accordance with the movement phases of an examined heart, in 
order thus to obtain tomograms of the heart in specific 
movement phases . 

[0024] — In accordance with the basic idoa at least one embodiment 
of the invention, the inventors also propose a CT unit for 
scanning an examination object that is provided with a detector 
array that is of planar design and has a multiplicity of 
distributed detector elements for detecting the rays of the 
beam from at least one focus, the at least one focus being 
movable relative to the examination object on at least one 
focal track that runs around the examination object and a 

detector array situated opposite, at least mo an a a 

device/method for collecting detector data, filtering and 3D 
back projection being provided, and ^feke — mcano a device /method 
for processing the measured data being fashioned in such a way 
that the method according to at least one embodiment of the 
invention as outlined above can be carried out. Said — The 
functional mcano aspects are — 

[0024] preferably implemented at least partially by programs or 
program modules. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0025] The invention is described in more detail below with 
reference to a n example preferred — exemplary — embodiment of a 
third-generation multirow CT unit suitable for carrying out the 
method according to at least one embodiment of the invention, 
and with the aid of the figures in which: 



Figure 1 shows multirow unit in a perspective illustration of 
the scanning unit in a schematic illustration of the 
peripheral equipment ; 

Figure 2 shows a longitudinal section through the multirow CT 
from Figure 1 ; 

Figure 3 shows illustrations of the recording geometry of a 
spiral CT unit; 

Figure 4 shows a schematic of the collection of measured data 
and of the approximate redundancy of measuring beams; 
and 

Figure 5 shows an illustration of the angular interval A(V) . 



DETAILED DESCRIPTION OF THE EXAMPLE EMBODIMENTS 

[0026] Figures 1 and 2 show a partially perspective 
illustration of a third-generation multirow CT unit suitable 
for carrying out the method according to at least one 
embodiment of the invention. The measuring arrangement (= 
Gantry) denoted by 1 has an x-ray source 2 with a beam 
diaphragm 3 in front thereof and near the source. The two- 
dimensionally constructed array with a number of rows and 
columns of detector elements forms the detector system 5 and is 
represented in section in figure 2 with a beam diaphragm 6 in 
front thereof and near the detector. 



[0027] For the purpose of better clarity, only eight rows Li to 
L Q of detector elements 4 are illustrated in figure 1. However, 
the detector system 5 can also have another, preferably 



8 



New PCT National Phase Application 



Docket No. 32860-000908/US 



greater, number of rows without departing from the scope of the 
invention. Another two-dimensional arrangement of the detectors 
is likewise also possible. 

[0028] The X-ray source 2 with the radiation diaphragm 3, on 
the one hand, and the detector system 5 with the radiation 
diaphragm 6, on the other hand, are arranged opposite each 
other on a rotary frame 7 in such a way that a pyramidal X-ray 
beam which, during the operation of the CT unit, emanates from 
the X-ray source 2 and is collimated by the settable radiation 
diaphragm 3 and whose edge rays are designated by 8, strikes 
the detector system 5. In the process, the radiation diaphragm 
6 can be set . to correspond to the cross section, set by means 
of the radiation diaphragm 3 , of the X-ray beam in such a way 
that, in accordance with different operating modes, only that 
area of the detector system 5 that is struck directly by the X- 
ray beam is exposed. In figures 1 and 2, only eight rows of 
detector elements 4 are used, the further rows indicated by 
dots are covered by the radiation diaphragm 6 and are therefore 
inactive . 

[0029] The X-ray beam has a cone angle (p which is the opening 
angle of the X-ray beam in a plane containing the system axis Z 
and the focus F. The opening angle of the X-ray beam in a plane 
lying at right angles to the system axis Z and containing the 
focus F (fan opening angle) is 2/?^. 

[0030] The rotary frame 7 can be set rotating about a system 
axis designated by Z by means of a drive device 22. The system 
axis Z runs parallel to the z-axis of a three-dimensional 
rectangular coordinate system illustrated in figure 1. 

[0031] The columns of the detector system 5 likewise run in the 
direction of the z-axis, while the rows run transversely with 
respect to the system axis Z. 
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[0032] In order to be able to bring an examination object, for 
example a patient, into the beam path of the X-ray beam, a 
bearing device 9 is provided, which can be displaced parallel 
to the system axis Z, that is to say in the direction of the z- 
axis, specifically in such a way that there is synchronization 
between the rotational movement of the rotary frame 7 and the 
translational movement of the bearing device, with the effect 
of producing a spiral movement of the focus and of a corotating 
detector relative to the examination object. The ratio between 
translational and rotational speeds can be set by a desired 
value for the feed H of the bearing device 9 being selected per 
rotation of the rotary frame. 

[0033] It is therefore possible for a volume of an examination 
object which is located on the bearing device 9, to be examined 
in the course of volume scanning, it being possible for the 
volume scanning to be undertaken in the form of spiral scanning 
in such a way that, with simultaneous rotation of the measuring 
unit 1 and translation of the bearing device 9, a multiplicity 
of projections from various projection directions are recorded 
by me an a — way of the measuring unit per revolution of the 
measuring unit 1. During the spiral scanning, the focus F of 
the X-ray source is moved relative to the bearing device 9 on a 
spiral path S about the examination object. 

[0034] M easured data read out in parallel from the detector 
elements of each active row of the detector system 5 during the 
spiral scanning, and corresponding to the individual 
projections in cone beam geometry, are subjected to 
digital/analog conversion in a data conditioning unit 10, are 
serialized and transmitted to an image computer 11. 

[0035] After the measured data have been preprocessed in a 
preprocessing unit 12 of the image computer 11, the resulting 



10 



New PCT National Phase Application 
Docket No. 32860-000908/US 



data stream passes to a volume image reconstruction unit 13, 
which uses the measured data to reconstruct images of the 
desired volume of the_— examination object by mcana — way of the 
computing operation described later. Both in the preprocessing 
unit 12 of the image computer 11 and in the volume image 
reconstruction unit 13, programs or program modules are stored 

by mcano — which the method according to at least one 

embodiment o f the invention can be carried out in order to 
generate CT images with the aid of a 3D image reconstruction 
method. 

[0036] The CT images are assembled from voxels (voxel = volume 
element) assembled like a grid, the voxels being assigned to 
the respective image volume, each voxel being assigned a CT 
number in Hounsf ield Units (HU) , and it being possible for the 
individual voxels to be displayed in accordance with a CT 
number/gray value scale in a gray value corresponding to the 
respective CT number. The image volume can be visualized in 
different ways. In a simple form, arbitrary planes can be 
displayed in an arbitrary orientation as tomograms. However, 
there are also more complex methods that visualize the entire 
volume. The shaded surface display (SSD) and volume rendering 

(VR) may be mentioned here by way of example. The SSD 
calculates for a settable threshold value an isosurface for 
which the image volume has the threshold value. The isosurface 
can be displayed on a screen with the aid of methods of 
computer graphics. In the case of VR, each voxel is assigned 
settable optical properties such as opacity and color in 
accordance with its value. Artificial views of the objects thus 
defined are calculated with the aid of methods of computer 
graphics . 

[0037] Images are calculated by a visualization unit 15 from 
the volume data reconstructed by the image reconstruction unit 
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13 and are displayed on a display unit 16, for example, a 
monitor connected to the image computer 11. 

[0038] The X-ray source 2, for example an X-ray tube, is 
supplied by a generator unit 17 with the requisite voltages and 
currents, for example the tube voltage U. 

[0039] In order to be able to set the latter to _ the 
respectively requisite values, the generator unit 17 is 
assigned a control unit 18 with a keyboard 19, which permits 
the values to be set as required. It is to be noted in this 
regard that in addition to the keyboard 19 it is also possible 
to provide other known input devices such as mouse, joystick or 
else voice input. 

[0040] In addition, the operation and control of the CT unit 
apart from this is carried out by moano use of the control unit 
18 and the keyboard 19, which is illustrated by the fact that 
the control unit 18 is connected to the image computer 11 . 

[0041] A mongst other things, the number of the active rows of 
detector elements 4, and therefore the position of the 
radiation diaphragms 3 and 6, can be set, for which purpose the 
control unit 18 is connected to the adjustment units 20 and 21 
assigned to the radiation diaphragms 3 and 6. In addition, the 
rotation time T can be set, which is the time needed by the 
rotary frame 7 for a complete revolution and which is 
illustrated by the fact that the drive unit 22 associated with 
the rotary frame 7 is connected to the control unit 18. 

[0042] The recording geometry of the spiral CT unit from 
figures 1 and 2 is illustrated schematically in figures 3 and 
4, the detector being reproduced only in a planar fashion here 
for the sake of simplicity. The radiation source moves with its 
focus F along a spiral track S of pitch h, which corresponds to 
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the bed feed, in a fashion rotating about an elongated 
examination object P, and in the direction of the system axis 
Z. In this process, the radiation intensity of the rays 
penetrating the examination object P is measured and collected 
on the detector D situated opposite. A volume V R to be 
reconstructed is embedded within the examination object P. 

[00 4 0] [0043] In order to reconstruct a volumetric image from 
the measured data, the latter must be filtered, and the 

filtered data must — subsequently be back-projected in three 

dimensions in order to generate a volumetric image of the 
examination object, the volumetric image then representing 
absorption values, obtained from the data, of the voxels V, 
belonging to the volume of the examination object, for the 
radiation of the beam. 

[0044] Thus, consideration of the reconstruction method shows 
that the cone beams from all focal positions that cut an image 
voxel V and are measured on the detector surface contribute to 
the reconstruction of the voxel V. The focal angular range A(V) 
with useful projection data is therefore a coherent angle 
interval, and is given by 

A(V)=[Oo(V) ,a x (V) ] . 

[0045] Figure 5 illustrates this angular interval A(V). It is 
assumed in this illustration without limitation of generality 
that the spiral track is traversed by focus and detector in a 
counterclockwise sense. The source position F{Oq) is 
distinguished in that the conical projection of the voxel V 
firstly strikes the detector surface D(a 0 ) corresponding to 
the source position F(cc 0 ), specifically the upper edge of the 
detector. In the further course of data acquisition, the 
conical projection of the voxel V migrates over the detector 
surface until it leaves the detector at the lower edge in the 
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case of the focal position F((Xi) . All the projections in the 
angular interval A(V) between (%> and Oi contribute to the 
reconstruction of the image voxel V. 

[0046] For example, the associated detector images can now be 
processed as follows: 

The first step is for the detector image to be subjected to a 
cosine weighting known per se, in order to take account of the 
oblique radiation present owing to the cone beam geometry. 

[0047] Subsequently, all the detector data are collected per 
voxel V considered, this being done along a straight line that 
runs through the cone beam projection of the voxel V and is 
aligned in the direction of projection t of the spiral tangent 
S t . If these data are not originally available, the required 
data are obtained from neighboring original data by mcano of via 
an interpolation method. 

[0048] This is followed by the approximate weighting for the 
purpose of compensating data redundancies. For this purpose, 
the projection angle 0 = a + P and the parallel coordinate p = 
R F . sin P (with R F being the radius of the focal track) can be 
determined from the focal angle a and the fan angle P, the cone 
angle being neglected. Two measuring beams S a and Sb are 
considered as redundant in this case when it holds that (0 a = 
2k-7i + 6 b and p a = Pt>) or (G a = (2k + l)-7t + 0 b and p a = -p b ) and k 
represents an arbitrary natural number. If redundant data are 
available, they can be multiplied by the generalized Parker 
weights- — sts — ars — described — ±& — E-8-K It is to be noted in this 
regard that the measuring beams defined as approximately 
redundant because of the different inclination to the z-axis 

(cone angle) are different and therefore not redundant. 
However, this approximation of the data redundancy can be 
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extended relatively far because of the generally relatively 
small cone angle. 

[0049] Ramp filtering with a suitable smoothing window can be 
used to filter the values, and the 3D back projection into the 
voxel V can be performed with the aid of a distance weighting 
known per se. 

[0050] It may be pointed out in addition that the method 
according to at least one embodiment of the invention is 

suitable both for CT units that are operated with 

simultaneously circulating focus/detector pairs . It is further 

suitable af*d — for CT units that have annularly arranged 

detectors that are fixed in the direction of rotation and i** 
the — cqog of w here ich only the_— focus rotates and the focus and 
detector can move in the z-direction. 
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[0052] Example embodiments being thus described, it will be 

obvious that the same may be varied in many ways. Such 

variations are not to be regarded as a departure from the 
spirit and scope of the present invention, and all such 
modifications as would be obvious to one skilled in the art are 
intended to be included within the scope of the following 
claims . 
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